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We consider in detail the quantum-mechanical problem associated with 
the motion of a one-dimensional particle under the action of the double-well 
potential. Our main tool will be the euclidean (imaginary time) version of the 
path-integral method. Once we perform the the Wick rotation, the euclidean 
equation of motion is the same as the usual one for the point particle in real 
time, except that the potential at issue is turned upside down. In doing so, 
our double-well potential becomes a two-humped potential. As required by 
the semiclassical approximation we may study the quadratic fluctuations over 
the instanton which represents in this context the localised finite-action solu- 
tions of the euclidean equation of motion. The determinants of the quadratic 
differential operators are evaluated by means of the zeta-function method. We 
write in closed form the eigenfunctions as well as the energy eigenvalues cor- 
responding to such operators by using the shape-invariance symmetry. The 
effect of the multi-instantons configurations is also included in this approach. 



* email : casahorr an@eupla. unizar . es 
tDFTUZ/00/09 



1 



I. INTRODUCTION. 



The semiclassical approximation to the euchdean path-integral represents one of the most 
valuable tools in quantum theory. By expanding around topological configurations one can 
obtain information about the tunneling phenomenon in physical systems ranging from non- 
relativistic one-particle quantum mechanics to relativistic gauge field theories. The basis 
of this approach relies on the so-called instanton which represents a localised finite-action 
solution of the euclidean equation of motion where the time variable is essentially imaginary. 
To sum up, we first find the appropiate classical configuration and subsequently evaluate 
the quadratic fiuctuations. The functional integration is solved in terms of the gaussian 
scheme except for the zero-modes which appear by virtue of the translational invariance 
of the system. It is customarily assumed the existence of collective coordinates so that as 
a matter of fact one performs the gaussian integration in the directions orthogonal to the 
zero-modes. 

To illustrate the usual instanton method in a simple model we resort to the double-well 
potential as first considered by Polyakov in his pioneering work on the subject. The 
famous problem of the level splittings is analyzed in terms of the instanton calculus which 
represents by itself an alternative to the WKB technique. From a physical point of view 
the result can be summarized as follows: the quantum mechanical amplitude consists on an 
exponential of the classical action of the topological configuration multiplied by a fiuctuation 
factor which takes into account the inverse square root of the determinant of harmonic 
eigenmodes. In this paper we will not try to expose new fundamental developments on the 
matter but rather focus the attention to the explicit evaluation of the determinants at issue. 

By analogy with the finite-dimensional case one can suspect that a functional determi- 
nant includes an infinite product of eigenvalues. As expected this expression will be highly 
divergent but in any case we may regularize it by means of the ratio of two determinants. 
When the eigenf unctions as well as the energy eigenvalues of the quadratic operators are 
known in closed form the zeta-function method serves to evaluate this quantity in a sys- 
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tematic way For the double- well potential it proves convenient to take advantage of 
the shape-invariance symmetry so that the only tools we need are supersymmetric quantum 
mechanics and the physics of the free-particle. 

In order to obtain an accurate description of the level splitting for a double-well potential 
(or double-humped once the potential is turned upside down) we introduce a chain of instan- 
tons and anti-instantons. If the distance between the centers of such configurations is large 
the interaction can be neglected. This framework represents the dilute gas approximation 
in quantum theory. The article is organised as follows. First of all we expose in detail the 
instantonic approach in non-relativistic one-dimensional particle mechanics. To illustrate 
the procedure in a simple but relevant context we consider the aforementioned double-well 
potential. We keep for the last section the analysis of the effect of the multi-instantons 
configurations. 



In this section we describe the instanton calculus for the one- dimensional spinless particle 
as can be found for instance in f^. For reference we assume that our particle moves under 
the effect of a confining potential V{x) so that the quantum model exhibits a pure discrete 
spectrum of energy eigenvalues. Whitout loss of generality we choose the origin of the energy 
such that the minimum (minima) of the potential satisfy V{x) = 0. Now taking the mass 
of the particle equal to unity, i.e. m = 1, the lagrangian which governs the behaviour of the 
model corresponds to 



If the particle is located at the initial time ti = — T/2 at the point Xi while one finds 
it when tf = T/2 at the point x/, the functional version of the non-relativistic quantum 
mechanics allows us to express the transition amplitude in terms of a sum over all paths 
joining the world points with coordinates (— T/2,Xj) and (T/2,Xf). The main physical 



II. INSTANTONS IN PARTICLE MECHANICS. 




(1) 
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ingredient consists in the way in which we weight the trajectories according to the imaginary 
exponential of the action. At this point it proves convenient to write the action S starting 
from the lagrangian L, i.e. 

fT/2 



S 



-T/2 



L(x, x) dt 



(2) 



so that the transition amphtude is 



< Xf\ exp{-iHT)\xi >= N{T) J [dx] expiS[x{t)] (3) 

It may be worth spelhng out that H represents the standard hamiltonian of the model 
while the symbol [dx] serves to indicate the integration over all functions which fulfill the 
aforementioned boundary conditions. The factor N{T) will be adjusted later on to make the 
whole integral finite and suitable normalised. As the hamiltonian H yields a pure discrete 
spectrum of energy eigenvalues, namely 



H\n >= En\n > 



(4) 



we have that 



< Xf \ exp{—iHT)\xi >= exp(— z-E^T) < Xf\n >< n\xi > 



(5) 



As this last expression is difficult to handle due to the oscillating character of the expo- 
nential, we perform the transition to the imaginary time. Making the substitution t — —it, 
known in the literature as the Wick rotation, one can take advantage of the euclidean for- 
malism where 
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iS[x{t)] 



T/2 
-T/2 



dr 



(6) 



Now it remains to identify the so-called euclidean action Se-, i.e. 
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5. 



T/2 
-T/2 



dr 



(7) 



The euclidean version of (§) corresponds to 



< Xf\ exp(-ifT)|x, >= N{T) J [dx] exp -5e[x(r)] (8) 

In order to explain the meaning of the semiclassical approximation we consider a function 
Xc(r) satisfying the boundary conditions at issue. Now we expand a general x{t) with the 
same boundary conditions according to 

x{t) = Xc{r) +Y,Cj Xjir) (9) 
j 

where Xji^r) constitute a complete set of orthonormal functions 

fT/2 

/ Xj{T) Xfc(r) dr = 5jk (10) 

vanishing at the boundary 

x,(±r/2) = (11) 

In addition the [dx] can be chosen as 

= n 4= (12) 

The semiclassical approximation takes for granted that xdr) represents a stationary 
point of the euclidean action, thus satisfying the equation 

^ - V'i^ (13) 

where the prime denotes as usual the derivative with respect to the spatial coordinate. We 
recognize in ([T3|) the euclidean equation of motion for our particle once the potential has 
been turned upside down. Now the crucial feature is the analysis of the second variation of 
the action with respect to Xc(t). (There is no linear variation due to the stationary point 
character of Xc{t)). In doing so we find that 

rT/2 

-T/2 



.T/2 

Se[Xc{T) + 6x{t)] = Seo+ / 

J -Til 



I d^ 1 



dr (14) 



being Seo the classical action associated with the configuration Xc{t). The procedure becomes 
much more accesible in terms of a complete set of eigenfunctions (eigenvalues) of the so-called 
stabihty equation 



^^v,{r) + V"[x^{T)]v,{T) = e,v,{T) (15) 



so that the family just anticipated in ([T0|) is obtained by means of ([15D . The problem 
becomes diagonal in this new scheme so that 

^e = ^eo + ^ES%- (16) 
^ j 

Taking into account the standard gaussian integral, i.e. 

/oo / Ott 

exp{-Xzy2) dz = (A > 0) (17) 

the euclidean transition amplitude reduces to 

< Xf\ exp{-HT)\x, >= N{T) exp(-5,„) Hej'^^ (18) 

j 

Almost from the very beginning of the subject the product of eigenvalues has been written 



as 



^1/2 



(19) 



a notation which obviously comes from the finite-dimensional case. It should be emphasized 
that by now we assume that all the eigenvalues ej are positive to avoid problems of definition 
for (|19D due to the existence of zero- modes (more on this later). The physical meaning of 
(p!8|) can be summarized as follows: the transition amplitude consists on an exponential of 
the classical action of the configuration Xc{t) at issue multiplied by a factor which takes into 
account the inverse square root of the determinant of a second order differential operator. 

To close this section we focus the interest on the explicit evaluation of the euclidean 
transition amplitude for the well-grounded problem of the harmonic oscillator. To be precise 
we consider boundary conditions such that = 0. The purpose of the calculation is 

twofold. First of all we find physical arguments to fix the factor N{T). On the other hand 
the harmonic oscillator acts as reference to deal with the ratio of determinants when we go 
to more complicated systems. Let us explore therefore the case where 
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V{x) 



X 



(20) 



so that V"{x = 0) = z/^. If the potential is turned upside down the euchdean equation of 
motion dictates that the only dynamical possibility starting from = at = — T/2 is 
Xc(t) = whenever we require also that = when = T/2. Otherwise the particle 
scapes towards plus or minus infinity. As the classical action associated with this trivial 
path vanishes the semiclassical approximation reduces itself to 



< X/ = 0| exp{-HhoT)\xi = >= N{T) jz^et 



dr^ 



-1/2 



(21) 



For the explicit evaluation of ( |21| ) we follow the pattern exposed in To start from 
scratch we write a complete set of eigenf unctions of the quadratic operator whose determi- 
nant we desire to compute. To be precise 



Xj = cos Ar 



Xj = sin Ar 



(22) 
(23) 



Once we add the boundary conditions (see (|TT])) the eigenvalues of energy ej which 
correspond to (|2^) and (^Sj) are of the form 



e, = ^ + ^^ (J = l,2,...) 



(24) 



Now the principal tool for extracting the physical information hidden in the product of 
eigenvalues is the factorization of the determinant, namely 



N{T) {Det 



d^ 



dr'^ 



-1/2 



N{T) 





-1/2 - 


'00 -2^2 

n 7^2 





n 1+ 



j2^2 



-1/2 



(25) 



As can be seen easily the two first factors represent the euclidean transition amplitude 
for the free-particle so that 

-1/2 



N{T) 



00 -2 2 



2^2 



=< Xf = 0\exp{-HoT)\xi = > 



(26) 
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where Ho corresponds to 



Ho 



(27) 



Now it suffices to resort to a set of plane waves labeled by |A; > to obtain that 

-1/2 

= g (28) 



oo -2 2 



2^2 



A^(r) 

Using again the formula (p!7| ) we find 

N{T) 



oo -2 2 



L ^2 



-1/2 



(29) 



The last factor of (125|) can be transformed into a more suitable expression according 

toi 

oo / ^2 \ 



n 1+ 



sinh TTZ 



(30) 



where in our case it happens that 



z/T 



TT 



(31) 



With all these partial results at hand the final form for the amplitude written in (^ 
should be 



< X/ = 0| exp{-HhoT)\xi = >= 



7r 



(2 sinh z/T) 



-1/2 



(32) 



For reference we point out the behaviour of (^) when considering the limit T oo 
since this final step is crucial for instance in the analysis of the double-well potential to get 
valuable information about the splitting. To sum up 



< Xf = 0\exp{-HhoT)\xi = > 



TT 



1/2 



exp(-z/T/2) 



l + -exp(-2z/T) + .. 



(33) 



Although it is a simple matter to check in ( p2D how this semiclassical approximation is 
exact for the harmonic oscillator we prefer to go to more relevant physical systems where 
the euclidean version of the path-integral represents probably the most elegant tool to study 
the tunneling phenomenon. 
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III. THE DOUBLE- WELL POTENTIAL. 



Once we have understood the main properties of the instanton calculus for the one- 
dimensional particle we proceed to examine the method in a concrete example. For such a 
purpose one needs a particle moving in a potential well which has a set of degenerate minima. 
No matter what the specific details of the potential, it is the case that the euchdean equation 
of motion exhibits topologically nontrivial solutions with finite-action. As anticipated in the 
previous section these solutions connect adjoining minima of the potential at issue and make 
the transition through a region classically forbidden. Althoug a periodic potential should 
be an excellent benchmark in this context we consider now the double-well potential V{x) 
given by 



From a classical point of view there are two minima located at a;_ = —1 and x+ = 1 
as expected when dealing with a potential which enjoys the discrete symmetry x — > —x. 
One could just as well assume that cu^ » 1 so that the barrier is high enough to decompose 
the system into a sum of two independent harmonic oscillators widely separated from each 
other. Then the lowest state of the model as a whole has a twofold degeneracy. Qualitatively 
speaking the particle lives in the right well or in the left one while executes small oscillations 
around the points a;_ = — 1 or = 1. For a finite barrier however, the particle in either 
of the two wells has a non- vanishing transition amplitude to tunnel into the other well so 
that the individual wave functions for left-hand and right-hand harmonic oscillators will 
mix. To sum up, the discrete symmetry x —>■ —x is not spontaneously broken since the 
expectation value of the coordinate x evaluated for the ground state is zero as expected for 
a wave function which makes an even superposition of wave functions centered at either of 
the two left and right wells. 

Now we can carry things further and discuss how the tunneling phenomenon emerges in 
the euclidean framework. The analysis is based on the existence of a topological configuration 




(34) 
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which connects the two minima located at x = ±1. The solution of the euclidean equation 
of motion which makes the connection is precisely the instanton according to the name first 
coined in the mid-seventies. Following the conventional wisdom the instanton interpolates 
between adjoining minima of the potential and lifts the degeneracy of the classical vacua. 
The instanton only appears once the potential is turned upside down. In doing so the two- 
well model transforms into a double-humped potential. At first glance it seems physically 
relevant to take into account the transition between the points x_ = — 1 and x+ = 1 
whenever we have previously managed to obtain the topological solution in closed form. 

Therefore we are interested in a trajectory which starts at ti = —T/2 from the point 
Xi = —1 and reaches xj = 1 when tf = T/2. It is customarily assumed that T ^ oo since 
the solution of the problem for finite T is much more complicated (more on this later). This 
movement, where the particle itself takes an infinite time to climb up (or drop from) the 
small piece next to the top of the two potential mountains, can be understood in terms of the 
conservation law of the euclidean energy once the double-well is just reversed. Simply put, 
we solve the problem by integration of a first-order differential equation instead of resorting 
to the conventional euclidean equation of motion written in (|T^. As the euclidean energy 



is equal to zero for the motion of our particle, the form of the topological solution is quite 
simple, namely 

Xc{t) = tanh ^ ^ (35) 

where the parameter indicates the point at which the instanton makes the jump. Of course 
the so-called antiinstanton is obtained by means of the transformation r — r and allows 
the connection between Xi = 1 and Xf = —1. Now we can compute the action associated 
with Xc{t) according to (0) so that 

Seo = f (36) 

It may be interesting at this point to remind that actually we need configurations for 
which X is equal to ±1 at large but finite values r = ±T/2. However the topological solution 
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written in (|35D reaches the points ±1 at infinite times. Fortunately the difference is so small 
that it can be ignored since at the end of the procedure we are interested precisely in the large 
T limit. As required by the semiclassical approximation we must calculate the quadratic 
fiuctuations correction once the classical contribution is known. The question we wish to 
address now is the analysis of the transition amplitude between Xi = —1 and xj = 1 taking 
as reference the harmonic oscillator derived from 



V"(X = ±1) = Lj' 



(37) 



Our description takes over 



<Xf = l\ exp{-HT)\xi = -1 >= N{T) i^Det 



-1/2 



Det [-{dydT^) + V"[x,{T)]]' 



-1/2 



exp(-5'eo) 



(38) 



[ Det[-{dydT^)+LU^] 

where we have multiplied and divided by the determinant for the harmonic oscillator pre- 
cisely. Incorporating the explicit form of Xc{t), where we fix by now Tc = 0, the stability 
equation reads 

3^2 



d^ , , 



Xj{T) 



e.jXj\T) 



(39) 



2cosh2(cjr/2)_ 

and corresponds to a Schrodinger's equation with Posch- Teller potential. Although the pro- 
blem can be studied in more conventional forms we prefer to do that by means of the shape- 
invariance symmetry 0. Ultimately the only tools we need are supersymmetric quantum 
mechanics and the free particle. (The details can be consulted in Appendix A). Among 
other things one finds a zero-mode Xoij) which in principle could jeopardize the evaluation 
of the determinant. However this eigenvalue eo = comes by no surprise since it refiects 
the translational invariance of the system. In other words, there is one direction in the 
functional space of the second variations which is incapable of changing the action. As a 
matter of fact one can discover the existence of a zero-mode starting from (|T^). It suffices 
an additional derivative with respect to r to find that 
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is just the solution of ( ]T5D with eg = 0. Notice the normahzation of Xo(r) which is due 
precisely to the zero euclidean energy condition for Xc(r). The way out of this apparent 
cul-de-sac is simple. The integration over Co (see (^) becomes equivalent to the integration 
over the center of the instanton Tc- To fix the jacobian of the transformation involved we 
take a first change such that 

Ax(r) = x„(r) Ac„ (41) 

According to the general expression written in (P) we find that under a shift Atc the 
effect should be 

Ax(r) = -/Se^ Xo(r)Ar, (42) 
Now the identification between (^) and (^2[) yields 



dco = \J SeodTc (43) 

where the minus sign disappears since what matters is the modulus of the jacobian at issue. 
In doing so we have that 

\ Det [-[£/dT^) + V"[x,{T)]] \-''' ^ 
\ Det[-{d^/dT^)+uj^] J 



\ Det'[-{d^/dr^) + V"W{T)]] \-"' [SZ, 

\ Det[-{d^/dT^)+uj''] ] \ 2n ' ^ ^ 

where Det' stands for the so-called reduced determinant once the zero-mode has been re- 
moved. 

Different methods can be found in the literature to find the ratio of determinants written 



in (0). One takes firstly the system enclosed in a box of lenght T, thus avoiding the subleties 
associated with the continuous spectrum At the end of the whole procedure it suffices 
the limit T ^ oo to achieve a meaningful result. However we prefer from the very beginning 
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to work in open space and evaluate the ratio of determinants by using the zeta-function 
method. (The interested reader can find in the Appendix B a sketch of this technique as 
tailored to our needs). Going back to (R9|) and (H), it suffices to make the change of variable 



LUT 



(45) 



to recognize the presence of the ratio of determinants (see Appendix B) 

Det' O2 



togheter with a global factor (5 given by 



Det P2 



(46) 



(47) 



As regards the non-zero spectrum of the operator O2 we have a discrete level with Ei = 3 
while the energy of the scattering states corresponds to 



Ek = k^ + A 



for eigenfunctions derived from the plane waves according to 



(48) 



Aliz) A\iz) 



exp{ikz) 



2tx 



VF+4 VW+i 

The explicit form of the operators (-2) and ^1(2;), namely 

d 



(49) 



dz 



+ tanh z 



d_ 

dz 



A\ = — ^ + 2 tanh z 



(50) 



(51) 



allows us to find that 



(t>2,k{z) 



cosh z 



— 3ik tanh z 



exp{ikz) 



(52) 



/27f vW+4: vF+T 

As a natural consequence of ( |5^ the regularized spectral density pr{k) of the problem 
reads 



13 



Pr{k) 



3(P + 2) 



(53) 



7r(P + 4)(P + 1) 

In this scheme the suitable zeta-function C,r{s) to evaluate the ratio of determinants of 
3) can be understood as 



This information is sufficient to write 

1 f°° r , ^ , 3 /•°° (A;2 + 2)exp[-(/c2 + 4)/i] 

^'^^'^ = W)Jo ' '^r^~''^~^L (P + 4)(F + 1) 

which, in turn, gives rise to ^ 

1 3 (A;2 + 2) 



(54) 



dk 



(55) 



Cris) 



dk 



(56) 



3^ ttJ-oo (P + 1)(P + 4)^+1 

Breaking the integral of (^) into more simple components it is not difficult to obtain 
Cr{s) in terms of Gamma and Hypergeometric Functions, namely 

1 3 1 r(s + |) 3 1 r(s + |) / 3 3\ 
^'^'^ " 3^ " 7^22^+1 r(s + l) " 7^22^+3 T{s + 2) ^ (^1' ^ + 2' " + 2' 4 j (^7) 

Taking into account that |^ 



Fll,-,2,- 
'2' '4 



(58) 



the evaluation of Cr(s) at s = yields 



Cr(0) 



(59) 



Let us come down to the concrete details which allow the evaluation of Cr(0)- Our 
approach provides a result of the form 



C;.(0) = - ln3 + 81n2 - i - ^F'{s + 1 1, s + 2, ^)|,=o 



(60) 



once one incorporates that 0] 



T\-) = -V^ (7 + 21n2) 



(61) 
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r'(i) = -7 



(62) 



r'(^) = -v^(| + ln2-l) (63) 



r'(2) = -7 + 1 (64) 

where 7 is the Euler's constant. In order to obtain the derivative of the Hypergeometric 
Function at issue we resort to the so-called integral representation, namely [|^ 

ns+li,s+2,^-) = J,^'J^^\, f\i-tr fi-^V'" dt (65) 



2' ' '4' r(i) r(s + 1) 7o ' 'V 4 

so that 

F'{s + ^,l,s + 2,^)U=o + + I, (66) 

for integrals of the form 

(67) 



/2 = " ln(l - t) dt (68) 



If we have that 



h - 3 (70) 



/2 = -ln- (71) 

3 3 ^ ^ 



/3 = -y + fln2 (72) 



then Cr(0) is obviously 
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C;(0) = 41n2 + ln3 (73) 

Under these conditions we can write finally the ratio of determinants R which appear in 
|), namely 

^ Det'[-{<f/dT^) + V"W{T)]] 

Det[-{(P/dT^)+uj^] ^ ' 

In a nutshell 

Combining now the results obtained along this section we write the transition amplitude 



anticipated in (|38D , i.e. 

< Xf = 1\ ex\;){—HT)\xi = —1 >= 



UJ 



(2sinha;T) ' y Seo \ — exp(— Sgo) cItc (76) 
71/ * V ^ 

Apart from the first factor, which represents the contribution of the harmonic oscillator 
of reference, we have a transition amplitude just depending on the point at which the 
instanton makes the jump. In accordance with the time interval T the result seems plausible 
whenever 

r- /6 

y^eo y- exp(-5'eo) T < 1 (77) 

a nonsense condition if T is large enough. But precisely in this regime we can accommodate 
configurations constructed of instantons and antiinstantons which mimic the behaviour of 
a trajectory just derived from the euclidean equation of motion. In doing so we get an 
additional bonus since the integration over the parameter is feasible. As a matter of 
fact the level splitting formula is deduced from the complete transition amplitude once we 
take into account the multi-instanton configurations. This serves to exhibit the instanton 
calculus as an alternative to the WKB approximation in this context. 
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IV. THE DILUTE-GAS APPROXIMATION. 



Although all the above calculations were carried out over a single instanton, it remains to 
identify the hypothetical contributions which incorporate the effect of a string of instantons 
and antiinstantons along the r axis. It is customarily assumed that these combinations of 
topological solutions represent no strong deviations of the trajectories just derived from the 
euclidean equation of motion without any kind of approximation. For background suppose 
we have j instantons and antiinstantons centered at points ri, r,- whenever 

T T 
- - < n < ... < r, < - (78) 

Being narrow enough the regions where the instantons (antiinstantons) make the jump, 
the action of the proposed configuration is almost extremal. We can carry things further 
and assume that the action of the string of instantons and antiinstantons is given by the 
sum of the j individual actions. This scheme is well-known in the literature where appears 
with the name of dilute gas approximation P|. 

It should be emphasized at this point that now we can compute amplitudes with closed 
paths with Xi = —1 = Xf for instance, so that the action at issue St will be an even multiple 
of the single instanton action, i.e. 

St ~ 2j S,, (79) 

As expected the complete transition amplitude between = — 1 and Xi = 1 incorporates 
the contribution 

St ^ (2j + 1) (80) 

In addition the translational degrees of freedom of the separated j instantons and anti- 
instantons lead to an integral of the form 

ujdTj \ ujdTj^i... / ujdTi = (81) 

-T/2 J-T/2 J-T/2 j\ 
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As regards the quadratic fluctuations around the j topological solutions we have now 
that the single ratio of determinants transforms into 



n 



Det [~{dydT^) +uj^ 



- exp(-c<;r/2) 



Det' [-{dydT^) + V"[xc{T)]] 



-1/2 



(82) 



Det[-(rfVrfr2) +cj2] J 

according to the limit of the factor associated with the harmonic oscillator when T is large. 
With all this information we can write the complete transition amplitudes for the double- well 
potential so that 



<Xf = l\exp{-HT)\x, = -!>=(- exp ( -cur/2) ^ ^ ^ 



TT 



(2j + l)! 



where d stands for the so-called instanton density given by 



6 / — 

d = ]l - y Seo exp(-S'eo) 



To sum up 



<Xf = l\exp{-HT)\xi = -!>= - exp(-cjT/2) sinh{ujTd) 



TX 



1/2 



Similarly 



<Xf = 1| exp(-i7T) |x. = 1 >= ( - ) exp(-.;T/2) 



so that 



< Xf = 1\ exp{—HT)\xi = 1 >= 



TT 



UJ 



IT 



j=Q 



1/2 



exp(— u;T/2) cosh.{uTd) 



(83) 



(84) 



(85) 



(86) 



(87) 



Taking the limit T — > oo in (pSD , we obtain the energy eigenvalues Eq and Ei of the first 
two levels of the double-well potential 



i?o ~ 2cijw — exp(— 2co'/3) 

2 V TT 



(88) 



Ei^ - + 2cjW-exp(-2cj/3) 
2 V vr 



(89) 



As anticipated in Section III the quantum mechanical tunneling transfers the wave func- 
tion from one well to the other, thus hfting the degeneracy of the classical vacua. To close, 
notice the way in which the energy eigenvalues depend on the barrier-penetration factor, i.e. 
the exponential of minus the classical action of the instanton. 
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V. APPENDIX A. 



In this appendix we include some results concerning the shape-invariance symmetry in 
one-dimensional supersymmetric quantum mechaniccs (susy qm for short). First of all one 
can assume the existence of two first-order differential operators A, which in terms of the 
superpotential function W{z) are of the form 

A^^ + W{z) (90) 
az 

At = - A + w{z) (91) 
az 

This enables us to write now a couple of hamiltonians and given by 

= A^A (92) 



To be precise 



H+ = AA^ (93) 



H^ = -^ + V^{z) (94) 



for 



V±{z) = W^{z)±W'{z) (95) 

where the prime denotes as usual the derivative with respect to the spatial coordinate. It is 
often the case to refer to the pair of potentials V_{z) and V+{z) as supersymmetric partners 
in accordance with the standard susy qm language. Now we look for the supersymmetric 
form of a standard quantum mechanical model associated with the hamiltonian Hy given 

by 

Hv^-^,+ V{z) (96) 
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where V{z) is some well-behaved function such that Hy has a non-empty discrete spectrum 
which is below a possible continuous spectrum. It seems plausible to identify Hy with the 
hamiltonian H of a standard susy qm model. Being e the energy of the ground state (f)o{z) 
of Hv we can write that 

Hv = H_ + e (97) 

so that e represents the role of the factorization energy while the superpotential W{z) which 
appears in the first differential operators A and is given by |^ 

In other words, we are dealing with a model where the ground state of Hy is nothing 
but the zero-mode of H_. Once written in the usual way, the two partner hamiltonians 
iJ_ and are essential isospectral: the ground state energy of vanishes and all other 
eigenvalues coincide with that of In such a case the hamiltonian Hy, defined as 

Hy = H+ + e (99) 

exhibits a spectrum which coincides with the set of energy eigenvalues of Hy except for e 
itself. It is obvious that by repetition of this method we can find a hierarchy of essential 
isospectral hamiltonians. 

At this point we can expose in brief how the so-called shape-invariance property allows 
to solve in a closed form the Schrodinger's equation. Let us assume the existence of a susy 
qm system described in terms of a superpotential W{z, Oq) which depends on the parameter 
tto- The partner potentials V±{z,ao) given by 

V±{z,ao) = W^{z,ao)±W'{z,ao) (100) 

are called shape-invariant if they are related by 

V+{z,ao) = V_{z,ai) + R{ai) (101) 
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where ai is a new parameter written in terms of Oq as ai = F{ao) while the residual con- 
tribution R{ai) is in fact independent of the variable z. To sum up, the shape-invariance 
property means that the partner potential V+{z,ao) can be considered as a new potential 
V-{z,ai) with superpotential W{z,ai) once the R{ai) has been subtracted. If the mapping 
ttg = F(as_i) may be iterated leading to a family of well-behaved superpotentials W{z,as), 
then the energy eigenvalues and their corresponding wave functions can be obtained in a 
simple way @]. 

As an example we focus our interest on the set of hamiltonians Oe given by 

^ dz"^ cosh^2; 

where i = 1, 2, although only the first members of the series are relevant in physics. It is 
the case that Oi can be factorized in terms of a superpotential W{z,i) like 

W{z,i) = itanhz (103) 

so that Oi = AIAi (notice that the energy of factorization e vanishes) for 

Ai = -^ + itanhz (104) 
dz 

AI = + itanhz (105) 
dz 

The shape-invariance condition appears once we write the partner hamiltonian Oe so 
that in this example the mapping between the old parameter i and the new one i reduces 
to 

i = i-l (106) 

The procedure may be iterated leading to a family of well-behaved superpotentials which 
allow us to extract the physical content of O^. The discrete spectrum includes a normalizable 
zero-energy mode 0£,o(-2) of the form 



J2(2i-1)\ 1 
- ^ (107) 
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which is below the set of states (f)e,m{z) given by 



^2(2£-2m-l)! 1 ^ ^ r 1 

'^^-(^^ = 2^-(£-m-l)! ynS^(E^ ^'^'^ ■ ■ ■ ^'"^^'^'^ 

with energies 

Ern^f-ie-mf 
ior m — 1, — 1. The continuous spectrum (f)e^k{z) corresponds to 

with energy eigenvalues 



Ek = k'^ + f 



27r 



and standard normalization as follows 



/oo 
4>i,k'{z) dz = S{k- k') 
-oo 
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VI. APPENDIX B. 



The first quantum correction to the classical term in the path-integral formalism is 
calculated by evaluating the determinant of a differential operator. As the determinant 
includes the product of the eigenvalues at issue the result would be in principle a terribly 
divergent expression. As far as we know Hawking 0] was the first to consider in a systematic 
way the zeta-function method as a useful technique for the study of these determinants. For 
reference, let us take a one- dimensional hamiltonian H with positive discrete eigenvalues Qj. 
If the corresponding eigenfunctions are (t)j{z) one can write that 

H(j)j{z) = aj(j)j{z) (113) 

Next we introduce the so-called zeta-function associated with H, i.e. 

Ch{s) = T.^ (114) 



Now if one notes that 

dCnis) 



^InOj exp(— slnoj) (115) 



ds ^ 
the determinant of H is given by 

Det H = exp -CniO) (116) 
This structure is enriched by the existence of the heat-kernel function G{z,w,fi), i.e. 

i 

which satisfies the heat difussion equation 

u ^ dG{z,w,n) 

H^G{z,w,iJ,) = (118) 

Of! 

with an initial condition like 



G{z, w, 0) = 6{z — w) 
24 



(119) 



On the other hand (nis) can be understood as the Melhn transformation of G{z,w,t) 
so that 

CHis) = -— /i^-irf/x / Giz,z,ij,)dz (120) 

1 [S) Jo J-oo 

To sum up, if we know the eigenf unctions (pjiz) and the eigenvalues aj of the hamiltonian 
H it suffices to insert such information in ( |117D to compute Ch{s)- Finally we use ( |llti| ) to 
obtain Det H. This zeta-function method makes it simple to study the scaling of the 
functional determinants we are involved with. Under a transformation of the form 

H ^ H = /3H (121) 

where jS represents a constant factor, one finds that 

Det H = 13^"^^^ Det H (122) 

as can be seen by inspection. However, this is not yet the end of the story since we need 
to adapt this method to the general case where the operator H also includes a continuous 
spectrum. To fix the ideas let us now sketch the way in which one can calculate the ratio of 
the determinants associated with Oi (see appendix A) and Pi which is given by 

p.--^.+e (123) 

Notice that Pi itself represents a free-particle comparison hamiltonian for Of,. To sum 
up we need to study 

where Det' Og denotes the reduced determinant obtained once the zero-mode has been 
explicitly removed. To handle the continuous spectrum we resort to the density matrix 
Eoi{k, z,w) written as [^] 

Eo,{k,z,w) = (j)lk{z) (j)e,k{w) (125) 

When going to P^ we get 
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exp [—ik{z — w)] 
2^ 



(126) 



as corresponds to a free-particle model. An additional subtraction prescription leads us to 
the regularized spectral density pr{k) expressed as 



Pr{k) 



[So, (k, z, z) - Tp, (/c, z)] dz 



(127) 



In this framework the appropiate zeta-function for the analysis of ( |124| ) would be 



Cr{s) = COf(s) - (p,{s) 



;i28) 



so that 



1 r°° f 



e-i 



- E [£2 - (£ - m)> 



m=l 



/oo 
prik) exp[-{k'^ + dk 
-oo 



(129) 
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